Magnetization of a half-quantum vortex in a spinor Bose-Einstein condensate 



Shoichi Hoshi and Hiroki Saito 
Department of Applied Physics and Chemistry, The University of Electro- Communications, Tokyo 182-8585, Japan 

(Dated: August 25, 2008) 

Magnetization dynamics of a half-quantum vortex in a spin-1 Bose-Einstein condensate with a 
ferromagnetic interaction are investigated by mean-field and Bogoliubov analyses. The transverse 
magnetization is shown to break the axisymmetry and form threefold domains. This phenomenon 
originates from the topological structure of the half-quantum vortex and spin conservation. 
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I. INTRODUCTION 

The spin degrees of freedom in an atomic Bose-Einstein 
condensate (BEC) allow for a variety of topological de- 
fects. Mcrmin-Ho [l| and Anderson- Toulouse Q vortices 
have been produced by the MIT group [3j using the topo- 
logical phase-imprinting method [4| . Recently, the Berke- 
ley group observed spontaneous formation of polar-core 
vortices in magnetization of a spin-1 87 Rb BEC [5|. In 
addition, several theoretical predictions have been made 
of various topological defects in spinor BECs, such as 
fractional vortices [1, 0, [1] , skyrmions @ , and knot struc- 
tures px| . 

Topological-defect structures in spinor BECs are 
closely related to the symmetry groups of spin states [III • 
For example, the symmetry group of the polar state 
(m = 0, where m denotes the magnetic sublevel of the 
spin) in a spin-1 BEC is U(l) x S 2 /Z2, and distinct con- 
figurations of the polar state can be specified by elements 
of this group. This symmetry group allows the topologi- 
cal structure of the half-quantum vortex [|| 0, [H • On the 
other hand, the symmetry group of the ferromagnetic 
state (|m| = 1) in a spin-1 BEC is SO(3). 

Let us consider a half-quantum vortex prepared in an 
antiferromagnetic BEC. We then consider a magnetic 
phase transition occurring from the antiferromagnetic to 
ferromagnetic phases. Since the symmetry group of the 
spin state changes from U(l) x S 2 /Z2 to SO(3) and the 
half-quantum vortex structure cannot exist for the lat- 
ter symmetry group, the ensuing magnetization dynam- 
ics are expected to break the symmetry and create non- 
trivial states. Such a change of spin texture associated 
with a change of the spin state symmetry group is the 
subject of the present paper. 

In the present paper, we study the magnetization dy- 
namics of the half-quantum vortex state produced in 
a spin-1 ferromagnetic BEC. The initial state can be 
prepared by, e.g., the phase-imprinting method using 
Laguerre-Gaussian beams [Hj]. We show that the half- 
quantum vortex develops into threefold magnetic do- 
mains through dynamical instability. This contrasts 
with the magnetization of a uniform polar state, in 
which a polar-core vortex or twofold domain structure 
is formed [13| . We study the dynamical instability using 
the Bogoliubov analysis, numerically in a 2D system and 



analytically in a ID ring. The threefold domain forma- 
tion can also be understood from the topological struc- 
ture of a half-quantum vortex and spin conservation, and 
its geometrical interpretation is provided. 

This paper is organized as follows. Section |TT] intro- 
duces the formulation of the problem and defines the no- 
tation used. Section HTT1 shows the magnetization dynam- 
ics of the half-quantum vortex and demonstrates three- 
fold domain formation. Section HVl details the Bogoliubov 
analysis to study the dynamical instability. Section [V] is 
devoted to the geometrical interpretation of the three- 
fold domain formation. Section IVT1 gives the conclusions 
to the study. 



II. FORMULATION OF THE PROBLEM 

We consider a BEC of spin-1 atoms with mass M 
confined in an axisymmetric harmonic potential V(r) = 
M[wjj_(x 2 +y 2 )-\-Lo 2 z 2 ]/2 that is independent of the mag- 
netic sublevels m of the spin. In the mean-field approx- 
imation, the condensate can be described by the macro- 
scopic wave functions ipmir) with m = —1, 0, 1 satisfying 



J dr ]T \^ m \ 2 = J drp - 



(1) 



where N is the number of atoms. The magnetization 
density is given by 



F ^ . ^Pmfrnva'^Prn'- 



(2) 



where / is the vector of the spin-1 matrices. The macro- 
scopic wave functions at zero temperature obey the three- 
component Gross-Pitaevskii (CP) equations, 
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where F± — F x ± iF y . The interaction coefficients in 
Eq. ((3]) are defined as 
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~M 3 : 



ci = 



AitTi 0,2 — do 
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with ao and 02 being the s-wave scattering lengths for col- 
liding channels with total spins and 2. For c\ < 0, the 
ferromagnetic state is energetically favorable, while for 
ci > the polar or antiferromagnetic state is favorable. 
In the present paper, we restrict ourselves to spin-1 87 Rb 
atoms, which have a positive cq and a negative c\ 



14] 



We assume that Tiuj z is much larger than other charac- 
teristic energies and the condensate has a tight pancake 
shape. The condensate wave function is therefore frozen 
in the ground state of the harmonic potential in the z di- 
rection and the system is effectively 2D. Integrating the 
GP energy functional with respect to z, we find that the 
2D wave function ip^ follows the GP equation having 
the same form as Eq. ([3]), where the interaction coeffi- 
cients Co and c\ are multiplied by [muj z / ^itTi)] 1 / 2 . We 
define a normalized wave function. 
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and normalized interaction coefficients, 
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with j ' = and 1. For example, using the scatter- 
ing lengths of a spin-1 87 Rb atom ao = 101. 8ae and 
a 2 = 100. 4ae [HI], where ae is the Bohr radius, and trap 
frequencies ujj_ = 2tt x 200 Hz and cu z — 2n x 4 kHz, the 
interaction coefficients become 

co ^ 0.167V, ci ~ -co/216. (7) 

We consider a half-quantum vortex state given by Q 



hqv 

hqv j 
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where r = (x 2 +y 2 ) 1 ^ 2 and 6 — a,rg(x+iy). The functions 
f±i(r) are stationary solutions of Eq. ^ satisfying 



r°° n 
J 2irr 2 \f ±1 {r)\ 2 dr = —. 



(9) 



From this condition, the state ([SJ) has an angular mo- 
mentum of Nh/2. Without loss of generality, we restrict 
ourselves to the upper sign in Eq. (|8|) unless otherwise 
stated. 

As several authors have discussed @, 0> @]> the half- 
quantum vortex has an interesting topological struc- 
ture. Equation ([5J is invariant under the transforma- 
tion exp(i</>/2) exp(if z (j>/2) exp(— 4>dg), where 4> is an ar- 
bitrary angle and f z — m for ip m . This indicates that 



spatial rotation by an angle cf> around the z axis is ac- 
companied by spin rotation by —<t>/2 with an additional 
phase factor exp(i(f>/2). Thus, for a rotation around the 
z axis by 2tt, the spin rotates only by — tt. 

Thermodynamic stability of the half-quantum vortex 
is studied in Ref. [16|. while a half-quantum vortex ring 
is discussed in Ref. [17|- Recently, it has been predicted 
that half-quantum vortices can be nucleated in rotating 
traps [3, [lj| and fluctuation-drive vortex fractionaliza- 
tion has been proposed in Ref. [2(| ■ 



III. MAGNETIZATION DYNAMICS OF A 
HALF-QUANTUM VORTEX 

In this section, we study the magnetization dynamics 
of the half-quantum vortex state ([5J) for a spin-1 87 Rb 
BEG 

The initial state is assumed to be the half-quantum 
vortex state ([SJ obtained by the imaginary-time propa- 
gation method. An experimental method to realize this 
initial state is discussed later. It follows from Eq. that 
when ipo is exactly zero as in Eq. ([8]) , "00 always vanishes 
in the mean-field evolution and no magnetization occurs 
even for a ferromagnetic interaction. We therefore add a 
small amount of initial noise to ipo, which triggers growth 
of the m = component. Physically, this initial noise 
corresponds to quantum and thermal fluctuations and 
experimental imperfections [H, [5l|, [52] • We set the noise 
as ipa = r i + i r 2 on each mesh point, where r\ and r-i 
are random numbers with uniform distribution between 
±10 -3 . For the imaginary- and real-time propagations, 
we employ the Crank- Nicolson m ethod with the size of 
each mesh being 0.05 ^Ti/ {muj±). 

Figure [T] (a) shows the time evolution of the density 
and phase profiles of each spin component and that of 
the transverse magnetization. At uo±t — 0, the trans- 
verse magnetization almost vanishes because of the initial 
state © with small noise added in the m = component. 
The m = component then grows in time and exhibits 
a threefold pattern, which leads to threefold magnetic 
domains as shown in Fig. [T](a) (u±t = 132). This three- 
fold domain formation is the main result of the present 
paper. We note that the magnetization in the three do- 
mains have different directions to cancel and conserve 
the total spin. Then, the population for each m oscil- 
lates as shown in Fig. Q] (b) due to the excess energy, and 
the appearance and disappearance of the threefold do- 
mains are repeated. At to±t ~ 300, the instability in the 
dipole mode becomes significant and the system under- 
goes a dipole deformation as shown in Fig.[T](a), which is 
followed by complicated dynamics. Throughout the dy- 
namics, the total density profile is almost unchanged and 
remains in the Thomas-Fermi distribution. The threefold 
domains are generated for both signs in the initial state 
([5]). We checked that no dynamics occur for an antifer- 
romagnetic interaction (ci > 0). 

We now discuss how to prepare the initial state in 
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FIG. 1: (color) (a) Time evolution of the normalized wave 
function ip m and transverse magnetization F+ — F x + iF y , 
where F x and F y are defined by Eq. (J2Jl with ip m replaced by 
tj} m . The size of each panel is 10 x 10 in units of y/Ti/(rmi)±). 
(b) Time evolution of the population in each component. The 
interaction coefficients are Co = 3000 and 5i = — So/216 (cor- 
responding to spin-1 87 Rb). 



an experiment. First we create a non-rotating BEC in 
the m = — 1 component, and then Gaussian and 

Laguerre-Gaussian beams propagating in the same di- 
rection arc applied. The frequencies of these beams are 
tuned to the Raman transition from the non-rotating 
m = — 1 state to the m = 1 state with a unit angular 
momentum [T^. We thus obtain the half-quantum vor- 
tex state given by 




^\e i9 sin(y4re- Sr " 


iP i ™ 1 cos(Are- Br ' 2 ) 



(10) 



where A and B are proportional to the intensity and 
width of the beams, respectively. For example, A — 
0.196(mw J _/n) 1 / 2 and B = O.OQ2mu±/n give density 
profiles similar to the initial state in Fig. [T] We have 
confirmed that the dynamics from this initial state are 
qualitatively the same as those in Fig. [TJ The pattern 
formation as shown in Fig. [T] (a) can be observed using 
the Stern-Gerlach separation and a nondestructive spin- 
sensitive measurement 1231. 



IV. BOGOLIUBOV ANALYSIS 
A. Numerical diagonalization 

The dynamics shown in Fig. [T] suggest that the half- 
quantum vortex state has dynamical instabilities. In this 
section, we perform the Bogoliubov analysis. 

We decompose the macroscopic wave functions ipm into 
the stationary state V'm v m Eq. © and small deviations 
Srp m from this state as 




e — i(/ii+|i_i)t/(2fi)^ o 
e -v-it/R(^ + ^_ 1 ) 



(11) 



where the chemical potential in each component is de- 
fined by 
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Here p and F z are given by Eqs. ((T|) and © with i/j m 
replaced by tp^ v . Substituting Eq. (JUJ) into Eq. ©, we 
obtain the Bogoliubov-de Gennes equations: 
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where we take only the first order of 5ip m . We note that 
both Eqs. (|13a[) and (|13b|) have a closed form within an 
angular-momentum subspace. Expanding dtp m as 



Slpr, 



^[afVH/^V)]^, (14) 



we find that otp only couples with (3[°} e , and oq' only 
couples with fl^lti a eZ\ \ and /S^Z/^- We therefore define 
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FIG. 2: Imaginary part of the Bogoliubov spectrum for 
the half-quantum vortex state of a spin-1 87 Rb BEC (5i = 
—co/216). The solid line is obtained by diagonalizing 
Eq. qi3ap . where the mode function has the form (|15p with 
1 = 2 and —1, corresponding to the threefold domain for- 
mation. The dashed line is obtained from Eq. ()13b[l . where 
the mode functions have the form (|16|1 with I — and 2, 
corresponding to the dipole deformation. 



the modes as 

8^ i =af\r)e ue e- i 
for the to — component and 
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for the to = ±1 components. 

We numerically diagonalize Eq. (fl"3)) using the method 
in Ref . [24[ . Figure [5] shows the imaginary part of the 
Bogoliubov spectrum for various values of cq and c\ with 
5i = —co/216. Diagonalizing Eq. (|13a|) for the to = 
component, we find that the excitation energies of the 
modes (|15[) with ^ = 2 and i = — 1 have an imaginary 
part for Co ^ 1450. Diagonalization of Eq. (|13b[) shows 
that the modes (p~6|) with £ = and I — 2 are dynamically 
unstable for Co > 0. The imaginary part for the to = 
modes exceeds that for the to = ±1 modes at Co — 
1870, and the growth in the to = component becomes 
dominant for cq > 1870. 

When the m = modes with 1 — 2 and £ = — 1 grow 
due to the dynamical instability, l^ol 2 becomes 
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We can show 



F + | 2 also has a similar form as a function of 9. 



Equation (jTTJ) indicates that these dynamically unstable 
modes generate threefold domains, in agreement with the 
result in Fig. [1] (a). The dynamically unstable modes 
of the m = ±1 components have angular momenta 1 ± 
1 in the to = 1 component and ± 1 in the to = — 1 
component. These modes therefore correspond to dipole 
deformation, which again explains the result in Fig.[T](a). 
Since the imaginary part of the m — excitation energy 
is larger than that of m = ±1 for cq = 3000, the threefold 
domains first emerge, followed by the dipole deformation. 



B. ID ring model 

For simplicity, we analyze a ID ring model in order to 
understand the dynamical instabilities in Fig. [2j 

We assume that the system is confined in a ID ring 
with radius R, and the effective interaction coefficients 
are denoted by Cq D and c\ D . For the stationary state 
V>„ qv in Eq. CP, we take 



V4 1 



hqv 



(18) 



where n = N/(2irR) is the atomic density and 9 is the 
azimuthal angle. The chemical potentials in Eq. (fT2"|) read 
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where 
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Substituting Eqs. (O and ((TSJ) into Eq. (fT5a|) gives 
.^ dSipo ( d 2 K 1D \ 1D i0 



Assuming that Sipg has the form 



5ip = a\'e 



(21) 



(22) 



we find that the coefficients o 
eigenvalue equations, 
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Diagonalizing these equations, we obtain the Bogoliubov 
eigenenergy for the to = excitation as 



huj 



(24) 
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For I = or 1, Eq. ((24)) is always real and there is no 
dynamical instability. For other values of £, the square 
root of Eq. (f24|) is imaginary when 2c\ D n < —K£(£ — 1), 
and the corresponding mode is dynamically unstable. 
The most unstable modes are £ — 2 and — 1. As in 
Eq. (|17|) . these modes correspond to the threefold domain 
formation. Thus, the transverse magnetization is most 
unstable against forming the threefold domains, which 



agrees with the 2D numerical result in Fig. [T] Similarly, 
the £ = 3 and —2 modes correspond to fivefold domains, 
the £ = 4 and —3 modes correspond to sevenfold domains, 
and so on. In general, dynamical instabilities forming j- 
fold domains with an odd integer j > 3 can exist. 

Performing the Bogoliubov analysis for 8ip±i in a sim- 
ilar manner, we obtain the eigenenergies as 
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where the corresponding mode has a form similar to 
Eq. (fT6|) with respect to £. For £ = 1, there is no dy- 
namical instability, since (^±1 have the same angular 
momenta as V±T in Eq. ([18]). For £ = and 2, Eq. ([25]) 
always has an imaginary part for cJ D < 0, in agreement 
with the 2D result in Fig. [2] The imaginary part is ex- 
panded as (cJ D c} D ) 1 / 2 n + 0(n 2 ). 



GEOMETRICAL MEANING OF THE 
THREEFOLD DOMAINS 




Now we consider the physical interpretation of the 
threefold domain formation. 
A spin state 



(Ci.Co,C-i) = (e ±l9 sin X ,0,cosx) 



has spin fluctuations as 



(26) 



FIG. 3: Lines indicating the directions of maximum fluctu- 
ations in the transverse magnetization for the half-quantum 
vortex state in Eq. (|26[) , where (a) and (b) correspond to the + 
and — signs, respectively. The dashed circles show schematics 
of the threefold domains and the arrows show the directions 
of magnetization in the domains. 
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sin 2 2 X , 

- [1 + cos(±6> + 2(j)) sin 2\] , 



(27) 
(28) 



where 9 is the azimuthal angle arg(ir + iy) and — 
f x cos 4> + f y sin (j>. The transverse fluctuation A/| then 
becomes maximum for <fi = ^9/2 and <fr = =F#/2 + 7r. The 
spatial distributions of the transverse fluctuation exhibit 
patterns, as shown in Fig. (3J where the direction of the 
line indicates the direction of the maximum transverse 
fluctuation. 

Spontaneous magnetization tends to occur in the direc- 
tions of large spin fluctuations, i.e., the directions of the 
lines in Fig. [3J For spontaneous magnetization, the total 
magnetization J drF must be conserved. From these two 



constraints, we understand the reason for the threefold 
domain formation. The arrows in Fig.[3]show examples of 
transverse magnetization satisfying the two constraints, 
where the magnetization in each domain occurs in the di- 
rection of the line and the sum of the three magnetization 
vectors vanishes. 

We note that continuous magnetization for all 9, as 
in the polar-core and Mermin-Ho vortices, is impossible, 
since the symmetry group of the ferromagnetic state is 
different from that of the spin state in Eq. (|26|) . The 
half-quantum vortex structure is peculiar to the latter 
symmetry group. Twofold domain formation is also im- 
possible, since spin directions at 9 and 9 + tt differ by 7r/2 
(not tt) and the total spin is not conserved. 

For j-fold domain formation, the center of each domain 
is located at 9 = 2np/j with p = 0, 1, • • • , j — 1. For each 
domain, there are two possible directions of magnetiza- 
tion, ^irp/j and ^fTrp/j+ir, where the =F signs correspond 



6 



to the ± signs in Eq. (|26|) . For the spin conservation to 
be satisfied, the sum of these magnetization vectors must 
vanish: 

i-i 

(e^ p/3 or e Tl(7rp/i+,r) ) = 0. (29) 

p=0 

For j = 1 and 2, Eq. (|29"|) cannot be satisfied. For j = 3, 
we find 1 + eT'W 3 ^)! e *27r/3 _ corresponds 
to the arrows in Fig. [3l In general, Eq. ([29]) can only be 
satisfied for an odd j > 3. This result agrees with that 
in Sec. ITVBl 

VI. CONCLUSIONS 

We have studied the dynamics of spontaneous mag- 
netization of a half-quantum vortex state in a spin- 
1 BEC with a ferromagnetic interaction. Solving the 
GP equation numerically, we found that the axisymme- 
try is spontaneously broken and the threefold magnetic 
domains are formed through the dynamical instability 
(Fig. [TJ. The critical strength of the interaction for 
the dynamical instability was obtained by the Bogoli- 
ubov analysis (Fig. [J). In order to understand the phe- 



nomenon in an analytic manner, we investigated the ID 
ring model and showed that the transverse magnetiza- 
tion is most unstable against forming the threefold do- 
mains among the j-fold domains with odd integers j ' > 3 
fSec. UVBj) . We provided a physical interpretation of the 
phenomenon based on the topological spin structure of 
the half-quantum vortex and spin conservation (Fig. [3|) . 

The half-quantum vortex in a spin-1 BEC is peculiar 
to the symmetry group that the spin state (|26[) pos- 
sesses. For the ferromagnetic interaction, in which the 
state ([26]) is unstable, the system exhibits nontrivial dy- 
namics, namely, threefold domain formation. We expect 
that various pattern formation phenomena may occur in 
magnetic phase transitions in spinor BECs containing 
topological structures, in which the symmetry groups of 
the spin states change in the phase transitions. 
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